Field Tuned Superconductor to Insulator Transitions in an Amorphous Film with an Imposed 

Multiply Connected Geometry 
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We have observed multiple magnetic field driven superconductor to insulator transitions (SIT) in amorphous 
Bi films perforated with a nano-honeycomb (NHC) array of holes. The period of the magneto-resistance, H = 
Hm = h/2eS where S is the area of a unit cell of holes, indicates the field driven transitions are boson 
dominated. The field-dependent resistance follows R(T) = Ro(H)e T °( H ' > ' T on both sides of the transition 
so that the evolution between these states is controlled by the vanishing of To — > 0. We compare our results 
to the thickness driven transition in NHC films and the field driven transitions in unpatterned Bi films, other 
materials, and Josephson junction arrays. Our results suggest a structural source for similar behavior found in 
some materials and that despite the clear bosonic nature of the SITs, quasiparticle degrees of freedom likely also 
play an important part in the evolution of the SIT. 
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Quasi two-dimensional electronic systems exhibit a wide 
range of phenomena including the quantum Hall effect 0], 
weak and strong localization 12[], and metal-insulator transi- 
tions Oil including the superconductor to insulator transition 
(SIT) 13El@]- These result from the enhancement of disorder, 
coulomb repulsion, and fluctuation effects on localization and 
superconductivity brought on by the reduced dimensionality. 
The aspects of these phenomena that are universal, or indepen- 
dent of the material comprising the systems, such as the exact 
quantization of the resistance plateaus in the quantum hall ef- 
fect [Q]] and the logarithmic decrease in the conductance of 
metal films with temperature Jst], have made them particularly 
intriguing. 

Despite many experiments, a universal description |01 for 
the SIT in amorphous films, where a superconductor is driven 
to an insulating state by increasing disorder or magnetic field, 
has not been established. While many experiments show scal- 
ing behavior around a critical point near the quantum of re- 
sistance for Cooper pairs, Rq = h/Ae 2 , qualitative differ- 
ences have emerged between the SITs of different systems. 
For instance, amorphous InOx and TiN films J! H EU show 
a spectacular evolution from a superconducting to an expo- 
nentially insulating state with increasing magnetic field. Be- 
yond this field tuned SIT, the magneto-resistance peaks at a 
value far in excess of both the normal state resistance, Rn, 
and Rq = h/4e 2 . Interestingly, not all amorphous materials 
show such dramatic behavior. In particular, amorphous MoGe 
dQ, Ta 0], Bi HQ, and Pb 0] show a resistance 
with a much weaker temperature dependence in their field in- 
duced insulating states and, at most, a very modest magneto- 
resistance peak iTTtIi . Some even exhibit a magnetic field in- 
duced metallic state El El El El El. While no complete 
explanation has taken root for this dichotomy, the most promi- 
nent explanations suggest that the spectacular behavior belies 
an underlying multiply connected geometry | T^, 2^1 ■ 



We have taken an amorphous film system (Bi) which ex- 
hibits the latter, weaker behavior and transformed it into a 
multiply connected system by patterning it with a regular 



nano honeycomb (NHC) array of holes. This transforma- 
tion results in a large magneto-resistance peak and a mag- 
netic field induced exponentially localized insulating phase, 
the salient features of the data on InOx and TiN. Moreover, 
these films exhibit multiple SIT's with increasing field imply- 
ing that Cooper pairs dominate both their superconducting and 
insulating ground states. 

These NHC films are produced within the ultra-high vac- 
uum environment of our dilution refrigerator cryostat. The 
templates used are aluminum oxide substrates specially an- 
odized to form the hole array 12111 . Those substrates used in 
this study had a center to center spacing a ~ 100 nm and a ra- 
dius Thole — 23 nm. AFM topographs show the substrate sur- 
face regularly undulates with a period ~ a, so that NHC films 
contain a regular (~ 20%) variation in the film thickness. Be- 
fore mounting in our cryostat these substrates are precoated 
with 10 nm of Ge and Au/Ge contact pads. Once at low tem- 
perature in the UHV environment of our dilution refrigerator, 
the substrates are coated with < 1 nm of Sb to ensure the sub- 
sequent Bi films' amorphous morphology. Through sequen- 
tial low-temperature evaporations of Bi we drive the system 
through a thickness tuned insulator to superconductor transi- 
tion without breaking vacuum or warming (see Fig. QJ). 

The four-point transport measurements were made using 
standard low-frequency ac techniques and confirmed through 
dc measurements of the IV characteristics with battery oper- 
ated electronics. The resistances of neighboring squares of the 
film were uniform to 5%. The magnetic field is supplied with 
a superconducting solenoid and is always perpendicular to the 
plane of the film. 

The main data presented come from the superconducting 
film closest to the critical point of the disorder tuned SIT of 
a series of films (see Fig. [TJi). It has a normal state sheet 
resistance, Rn, of 21 kfi and a mean field transition temper- 
ature, T c0 ~ 1.5 K JH. Below T c0 , the resistance, R(T), 
drops exponentially with an activated form. The thinner in- 
sulating films also exhibit activated resistances with an acti- 
vation energy of the opposite sign. Thus the superconducting 
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Figure 1: a) Resistance versus inverse temperature for the thickness 
tuned superconductor to insulator transition. Dashed lines are fits 
to R(T) = Roe T °/ T . To from most insulating to superconducting 
are 3.1, 2.65, 2.05, 0.37, and -1.76 K. b). Inverse temperature plot 
of the field driven 1ST bearing a strong resemblance to (a). Data 
in fields between H — and H = Hm/2 are nearly temperature 
independent. All dashed lines are fits to the same form with Tbs 
from top to bottom of 0.194, 0.108, 0.07, unfit, -0.04, -0.107, and 
-1.76 K. 



state emerges from exponentially localized state and the tran- 
sition can be characterized by a continuously varying energy 
scale 0. 

The initial magnetic field driven SIT appears in Fig. QJ>. As 
the field is increased, the low temperature dR/dT changes 
from superconducting-like, dR/dT > in H = 0, to 
insulating-like, dR/dT < for H > 0.15H M - The low tem- 
perature resistance continues to increase up to Hm/2. Here, 
Hm = h/2eS is called the matching field and S is the area 
of a unit cell of holes. Figure QJ shows the resistance can be 
parametrized by R(T) = R Q {H)e Toi - H ^ T at the lowest tem- 
peratures in all fields except for the that in H = 0.15-Hm- 
The energy scale, To, extracted from some of the fits shown is 
small compared to the ultimate temperature and so cannot be 
associated with a simple activated process. Nonetheless, the 
superconducting and most insulating (Hm/2) data are well 
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Figure 2: Activation energies as a function of magnetic field which 
show the 1ST is governed by a single parameter To — > 0. The box 
at To(H) — represents the data taken in H — 0.15Hm which is 
temperature independent. The dashed line is a guide to the eye. Inset 
shows Ro(H) depends approximately logarithmically on field across 
the SIT. The dashed line in the inset is a fit to this functional form. 



described by this form and it seems the most reasonable pro- 
cedure for interpolating between them. 

The energy scale characterizing the transition, To, varies 
smoothly through the SIT and changes sign for H ~ 0.15-Ha/ 
(see Fig. Thus it appears To — > sets the critical field, 
H c , while the value of the critical resistance is set by the other 
fit parameter, Ro(H = .15Hm) ~ 3 kfl. The prefactor Rq 
varies roughly logarithmically through the SIT (see Fig. [2] 
inset). 

Beyond Hm/2 dR/dT and Tq changes sign multiple times 
suggesting multiple field tuned SITs. Figure [3^ shows that 
dR/dT > at the matching fields 0, H M , and 2H M and 
dR/dT < at Hm/2, and '6Hm/2 as appropriate for su- 
perconducting and insulating behaviors, respectively. The 
magneto-resistance isotherms in Fig. [2^ further illustrate this 
multiplicity of transitions and the underlying matching field 
period. The inset of Fig. [2J5 provides the detailed shape of 
the first peak in the magneto-resistance. It appears symmetric 
about Hm/2. It is also apparent in Fig. [3^ that the magneto- 
resistance oscillations decrease in amplitude with increasing 
field. The "superconducting" state at Hm is more resistive 
than at H = and the insulating state at 3Hm/2 is less resis- 
tive than at Hm/2. These temperature dependences weaken 
further for subsequent SITs. In a field of a few Tesla, the asso- 
ciated magneto-resistance oscillations damp out completely. 

Beyond the oscillations, the films become increasingly in- 
sulating with magnetic field. Figure [4] also shows that the 
R(T) in higher fields maintain an activated form (at least up 
to H = 4 T) and the reentrance persists. Remarkably, the 
fit prefactor remains unchanged in a 4 T field from the value 
attained in Hm /2, nearly 20 times lower in field. 

The periodic appearance of "superconducting" states at 
matching fields set by the charge 2e flux quantum implies that 
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Figure 3: a) R(T) for half-matching and matching fields. Those data 
in matching fields show the tendency to superconduct as T — > 
showing there are at least three transitions between superconducting 
(dR/dT > 0) and insulating (dR/dT < 0) behavior with increasing 
field in the NHC film. Dashed lines are fits to R(T) = R e {T ° /T) 
with T from top to bottom of 0.194, 0.157, -0.03, -0.138, and -1.76 
K b) Magneto-resistance isotherms at a variety of temperatures. Inset 
shows a continuous field sweep over one period on a linear resistance 
scale. 



Cooper pairs are the primary charge carriers throughout the 
field driven SIT. Thus, bosonic degrees of freedom dominate 
this transition. Moreover, comparison of the disorder and field 
tuned transitions suggest that similar bosonic physics operates 
in the two cases. There are three common features. Each ex- 
hibits activated resistances on both sides of the SIT and quasi- 
reentrant minima in the R(T). Moreover, each can be char- 
acterized by a vanishing energy scale, To, at the critical point 
implying that the system is metallic there. These similarities 
suggest these SITs involve common mechanisms and that the 
magnetic field induces fluctuations of the same character as 
those induced by decreasing the thickness of NHC films. 

The disorder and field tuned SITs of similar unpatterned 
films dill differ qualitatively from the SITs of NHC films 
in a manner suggesting that fermionic, rather than bosonic de- 
grees of freedom dominate their behavior. The R(T) of unpat- 
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Figure 4: Resistance behavior in high magnetic field for the H — 
superconducting film. At high fields the resistance is activated and 
the magnetoresistance positive suggesting the same transport mech- 
anisms even when oscillations are small. Tb's from top to bottom are 
0.69, 0.51,0.194,-1.76 K. 



terned insulating films increase logarithmically rather than ex- 
ponentially with decreasing temperature and they exhibit lit- 
tle, if any, reentrant behavior. Moreover, H c for unpatterned 
films approaches the expected upper critical field, H c2 , or the 
field at which the amplitude of the order parameter goes to 
zero. By contrast, H c <C H c2 for NHC films, which al- 
lows multiple field tuned SITs. The above unpatterned film 
characteristics suggest that single electron motion dominates 
their insulating transport and order parameter amplitude fluc- 
tuations and pair breaking are strong in their superconduct- 
ing phase. Thus, fermionic degrees of freedom influence both 
sides of their SIT. 

On the other hand, the NHC films SITs qualitatively re- 
semble those of Josephson Junction arrays (JJA). These litho- 
graphically prepared systems also show multiple SITs and 2e 
periodic magneto-resistance oscillations provided the ratio of 
the Josephson coupling, Ej, to the charging energy, Ec, is 
nearly critical (i.e. Ej/E c ^ 1) JJAs also exhibit 



minima in their R(T) and activated resistances 124 12611 as 
do NHC films. While SITs do not occur at fractional values 
of H/Hm for NHC films as they do in JJAs, it is likely that 
quenched disorder in the NHC lattice smears these out. The 
SITs of JJAs are usually discussed in terms of the bose con- 
densation of vortices J3]. This picture seems to account for the 
R{T) over a limited temperature range JU It fails, however, 
for the lowest temperature data as the R(T) become temper- 
ature independent in the T — > limit 12411 . This effect seems 
to require coupling to quasiparticle degrees of freedom in the 
array whether intentionally introduced through shunt resistors 
lH or not J3. 

While the similarities above suggest that JJA physics drives 
the field driven SIT in the NHC films, there are some chal- 
lenges to a direct comparison. First, the NHC films do not 
contain high quality tunnel junctions. All portions of NHC 
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films are continuous, strongly disordered films without any 
high quality oxides. Second, the estimated charging energy 
using the nodes as the superconducting islands is far too large 
to compete with Ej ITM . Third, the R(T) do not scale well 
around the critical point of the SIT as they do for the SITs 
of JJAs. We suggest that a model including quasi-particle de- 
grees of freedom in addition to Bose degrees of freedom may 
be more appropriate for the NHC films ll20ll . Virtual quasi- 
particle processes can screen Coulomb interactions to renor- 
malize the charging energy to a value lower than the geometry 
implies 125l 12711 . Moreover, quasi-particle motion can provide 
a dissipative channel that could lead to metallic behavior and 
deviations from scaling II ill . 

The discussion above strongly suggests that patterning 
amorphous films into a multiply connected geometry pro- 
motes the influence of bosonic degrees of freedom. Recently, 
in fact, two models invoking such a geometry have emerged to 
explain the activated resistances and giant magneto-resistance 
peak near the SIT of InOx and TiN 18, 1(1 28] films. In one 
picture l2(ill . Cooper pair transport between islands 1E9I [30h 
takes place in parallel with quasiparticle transport in the nor- 
mal portions of the film. For H < H pea k, the magnetic field 
reduces phase coherence between islands, increasing the re- 
sistance. Quasiparticle transport is dominant for H > H pea k 
as the field decreases the number and size of islands so that a 
peak is produced when the two mechanisms balance. In the 
second 11911 . the film is treated as an array of SQUIDs. Here 
the peak is the result of the magnetic field tuning Ej — * 0. 
The position of the peak is then determined by the effective 
area of a squid loop in much the same way as the position of 
the first peak in Fig. is determined by the area of a unit 
cell of holes. Applying these ideas to NHC films seems to re- 
quire elements of both explanations. The position and period- 
icity of the magnetoresistance peak is best explained through 
Ej(H) oc cos(H) while the nearly metallic behavior of the 
R(T, H) (0 < H < H M /2, H M , and 2H M ) seems to require 
dissipation. In each case, the inherent order parameter mod- 
ulation fl9l I20I1 is probably enhanced in NHC films through 
their patterning and small thickness variations. NHC film be- 
havior, therefore, seems to provide a bridge between that of 
MoGe 10,0 and Ta on the one hand and TiN 0] and 
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InOx II10H on the other. 

We have measured multiple field driven SITs in an amor- 
phous Bi film patterned with a nano-honeycomb array of 
holes. Each SIT is boson dominated, however, a description 
in terms of bose condensation falls short. Here, a full explana- 
tion would seem to require Josephson physics and quasiparti- 
cle degrees of freedom in an effectively islanded system. The 
activated resistances ifioll and the magneto-resistance peak 
ifsLIToh that NHC films share with TiN and InOx and the nearly 
metallic behavior shared with MoGe | IT , 12 ] and Ta J3l in- 



dicate that NHC films occupy a space between these two ex- 
tremes. Moreover, it seems a dissipative coupling and the type 
and degree of islanding interpolates between them. 
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